AN OVERDETERMINED PROBLEM IN POTENTIAL THEORY 
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Abstract. We investigate a problem posed by L. Hauswirth, F. Helein, and F. 
Pacard 1191 , namely, to characterize all the domains in the plane that admit a 
"roof function", i.e., a positive harmonic function which solves simultaneously 
a Dirichlet problem with null boundary data, and a Neumann problem with 
constant boundary data. As they suggested, we show, under some a priori 
assumptions, that there are only three exceptional domains: the exterior of 
a disk, a halfplane, and a nontrivial example found in 1 191 that is the image 
of the strip [9£| < tt/2 under ( -5- ( + sinh(C). We show that in R 4 this 
example does not have any axially symmetric analog containing its own axis 
of symmetry. 



1. Introduction 

In |19| . the authors have posed the following problem: find a smooth bounded 
domain f2 in a Riemannian manifold M g with metric g, such that the first eigen- 
value Ai of the Laplace-Beltrami operator on ft has a corresponding real, positive 
eigenfunction u\ satisfying ui — 0, = 1 on the boundary of Q. Any such domain 
is called extremal because it provides a local minimum for the first eigenvalue Ai 
of the Laplace-Beltrami operator, under the constraint of fixed total volume of f2 
(see [H] and references therein). 

In special cases one can find a sequence of extremal domains {Qt} with increasing 
volumes, such that the limit domain f2 = f2t_>oo is unbounded, and its first eigen- 
value vanishes as t — > oo. This limit extremal domain is then called exceptional, and 
the corresponding limit function (tti t)t-><» — > it is a positive, harmonic function on 
n which solves simultaneously the overdetermined boundary value problem with 
null Dirichlet data and constant Neumann data. 

The problem of finding exceptional domains in WL n and their corresponding func- 
tions u (called "roof functions by the authors of 19J ) is a nontrivial problem of 
potential theory. There is no obvious variational principle to use, on the one hand 
because is unbounded (so the Dirichlet energy of u [U Ch. 1] will diverge), and, on 
the other hand, because the constant Neumann data constraint is not conformally 
invariant. 

In the absence of a suitable variational formulation, we may interpret the scaling 
t — > oo described above as a dynamical process, in which the pair (Tt tl u t ) evolves 
so that the limit t — > oo solves the overdetermined problem. In other words, we can 
turn this observation into a constructive method for finding (building) exceptional 
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domains. In order to do this, it is helpful to note that, upon compactification of the 
boundary dfl (with metric da 2 ), the pair (fi, u) with flat metric becomes conformal 
to the half-cylinder AT := R+ x dfl, with metric 

ds 2 = e~ 2u {du 2 + da 2 ). 

Under this reformulation, scaling of (fit, Ut)t-Hx> becomes equivalent to scaling of the 
metric structure given above, defined over the fixed space Af . This is reminiscent of 
the Ricci flow, in which the metric structure g evolves with respect to a deformation 
parameter feM according to the equation 

d 9ij _ OR 
It ~ ■ h 

with the right side of the equation given by the covariant Ricci tensor. It is known 
[37] that for the case of a two-dimensional manifold, with metric given by ds 2 = 
e~ 2u (dx 2 + dy 2 ), the Ricci flow equations reduce to a single nonlinear equation 

du ~ 

(since in two dimensions the Riemann tensor has only one independent component). 
This is a heat equation with the generator given by the Laplace-Beltrami operator 
corresponding to the metric ds 2 . Therefore, if there is a stationary solution — > 
as t — > oo, it will correspond to the scaling of the first eigenvalue Xi(t) — > and, 
by conformally mapping back Af using the solution u(t — > oo), we will obtain the 
solution (ft, u). 

In other words, we can summarize this constructive method for finding excep- 
tional domains in M 2 as follows: starting from a 2-dimensional Riemannian manifold 
with finite volume and metric encoded through the positive real function u, and 
boundary set defined via u = 0, consider the time evolution given by the Ricci flow, 
without volume renormalization. Then |37| the manifold will remain Riemannian 
at all times, and in the t — > oo limit the function u will become a solution of the 
nonlinear Laplace-Beltrami equation. Furthermore, if u remains finite everywhere 
in the domain, then it is harmonic and satisfies both Dirichlet and Neumann con- 
ditions at all finite boundary components, so it is a solution for the overdetermined 
potential problem. Considered together with the (boundary) point at infinity, the 
manifold is equivalent [3D] to a pseudosphere (flat everywhere except at the infinity 
point, with overall positive curvature). (We wish to emphasize that there is no 
reason to assume that such constructive methods would be exhaustive.) 

Thus, so motivated, it is natural to try to characterize exceptional domains in flat 
Euclidean spaces. The authors in [TH] suggested that in two dimensions there are 
only three examples: a complement of a disk, a halfplane, and a nontrivial example 
obtained as the image of the strip |3f | < 7r/2 under the mapping ( — > ( + sinh(£). 
They posed as an open problem to determine if these are the only examples [191 
Section 7] . (They gave some evidence by characterizing the halfplane under a global 
assumption on the gradient of the roof function (XSJ Prop. 6.1].) They also posed 
the problem of finding nontrivial examples in higher dimensions and suggested the 
possibility of axially symmetric examples similar to the nontrivial example in the 
plane [19 Remark 2.1]. 

We address both of these problems. The paper is organized as follows. In Section 
2, we review the theory of Hardy spaces in order to address a subtlety that arises 
in connection with the regularity of the boundary of an exceptional domain. This 
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leads us to assume in our theorems that the domain Q is Smirnov. In Section 3, we 
characterize exteriors of disks as being the only exceptional domain whose boundary 
is compact. In Section 4, we establish a connection between the "roof function" 
of an exceptional domain and the so-called Schwarz function of its boundary, and 
we also show that the boundary of a simply connected exceptional domain Q can 
pass either (i) once or (ii) twice through infinity. In Section 5, we show that Case 
(i) implies that Q is a halfplane. In Section 6, we show that Case (ii) implies that 
f2 is the nontrivial example found in [191 Section 2]. In each of these theorems we 
assume that D, is Smirnov, but we allow the roof function to be a weak solution 
merely satisfying the boundary conditions almost everywhere. 

In Section 7, we extend the result of Section 3 to higher dimensions. In Section 
8, we show that the nontrivial example from Section [6] does not allow an extension 
to axially symmetric domains in four dimensions, contrary to what was suggested 
in [ini Remark 2.1] (and we conjecture that this example has no analogues in any 
number of dimensions greater than two). 

Sections 3 through 6 together partially confirm what was suggested in [191 Section 
7] under some assumptions on the topology of f2 and assuming that Q is Smirnov. In 
Section 9, we give concluding remarks including a conjecture that, up to similarity, 
there are only three finite genus exceptional domains. The additional assumption of 
finite genus is due to a remarkable example of an infinite-genus exceptional domains 
that appeared in the fluid dynamics literature [4]. See Section 9 for discussion. 

Remark 1. After this paper was submitted, Martin Traizet announced a more 
complete classification of exceptional domains [33] after developing an exciting new 
connection to minimal surfaces. He characterized the three examples among those 
having finitely many boundary components. M. Traizet's preprint that appeared 
while we have been revising the previous version of our paper finds a new beautiful 
connection of the problem with the theory of minimal surfaces. From this point of 
view, he noticed the above-mentioned family of infinite genus examples [4] and also 
characterized them among periodic domains for which the quotient by the period 
has finitely many boundary components [381 Theorem 13]. For this latter result, he 
invokes a powerful theorem of W. H. Mccks and M. Wolf. Our methods mostly rely 
on classical function theory (H p spaces) and potential theory and in most parts are 
different from Traizet's. Interestingly, as Traizet notes in his preprint [38], Remark 
5], if one could prove his Theorem 13 by only invoking pure function theory this 
would give (via Traizet's results) a new and independent proof of the Meeks-Wolf 
result from minimal surfaces. An attractive challenge! 

Acknowledgement: The authors are indebted to Dimiter Vassilev for bringing the 
article [19j to their attention. We wish to thank Alexandre Eremenko for sharing an 
improved proof of Theorem |4.2| and Arshak Petrosyan and Koushik Ramachandran 
for pointing out the example of a cone as an exceptional domain. We also wish 
to thank Martin Traizet for helpful discussion regarding his preprint. The two 
first named authors acknowledge partial support from the NSF under the grant 
DMS-0855597. 

2. Classical vs. Weak Solutions, Regularity of the Boundary, and 

Hardy Spaces 

From the rigidity of the Cauchy problem, one might expect to obtain, "for 
free", regularity of the boundary of an exceptional domain (as is often the case 
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for solutions of free boundary problems). Unfortunately, the problem at hand is 
complicated by a remarkable family of examples with rectifiable but non-smooth 
boundaries, a.k.a. non-Smirnov domains - cf. [121 Ch. 10]. This results in adding a 
Smirnov condition to the assumptions on the domains if we desire to consider "weak 
solutions", i.e., harmonic "roof functions" satisfying the Dirichlet and Neumann 
boundary conditions almost everywhere with respect to the Lebesgue measure. 

In order to address this subtlety, we first give some background from H p theory 
- cf. [12 for details. 

An analytic function / : D — > C is said to belong to the Hardy class H p , < p < 
oo, if the integrals: 

^ \f(re i8 )\Pde 

o 

remain bounded as r — > 1. 

Recall that a Blaschke product is a function of the form 



B(z) = z™TT— - 

- LJ - On 1 



a n z 



where m is a nonnegative integer and — \a>n\) < 00 • The latter condition 
ensures convergence of the product (See Theorem 2.4 in [12]). 

A function analytic in D is called an inner function if its modulus is bounded by 
1 and its modulus has radial limit 1 almost everywhere on the boundary. If S(z) is 
an inner function without zeros, then S(z) is called a singular inner function. 

An outer function for the class H p is a function of the form 

f 1 f 27r e lt + z 1 
(2.1) F(*) = e^exp|— J _±_i 0g ^(i)dtj , 

where 7 is a real number, ip(t) > 0, \ogip(t) £ L 1 , and ip(t) £ L p . 

The following theorem [12, Ch. 2, Ch. 5] (also cf. [TC]) provides the parametriza- 
tion of functions in Hardy classes by their zero sets, associated singular measures, 
and moduli of their boundary values. 

Theorem 2.1. Every function f(z) of class H p (p > 0) has a unique (up to a 
unimodular constant factor) factorization of the form f(z) = B(z)S(z)F{z) , where 
B{z) is a Blaschke product, S(z) is a singular inner function, and F{z) is an outer 
function for the class H p . 

Suppose Q is a Jordan domain with rectifiable boundary and / : D — > Q is a 



conformal map. Then /' £ H 1 by Theorem 3.12 in [TJ]. By Theorem 2.1 /' has 
a canonical factorization f'{z) = B(z)S(z)F(z) 1 and since / is a conformal map 
/' does not vanish, so f'(z) — S(z)F(z). Then £1 is called a Smirnov domain 
if S(z) = 1 so that f'(z) = F(z) is purely an outer function. This definition is 
independent of the choice of conformal map. 

There are examples of non-Smirnov domains with, as above, f'(z) = S(z)F(z), 
but now F{z) = 1 and the singular inner function S(z) is not constant. Such 
examples were first constructed by M. Keldysh and M. Lavrentiev [53] using com- 
plicated geometric arguments. Their existence was somewhat demystified by an 
analytic proof provided by P. Duren, H. S. Shapiro, and A. L. Shields [13]. Like the 
disk, such a domain has harmonic measure at zero (assuming /(0) = 0) proportional 
to arc-length. Thus, its boundary is sometimes called a "pseudocircle" . 
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Similarly, there are "exterior pseudocircles" , arising as the boundary of an un- 
bounded non-Smirnov domain [2 2) for which the harmonic measure at infinity is 
proportional to arclength, and thus Green's function with singularity at infinity 
provides a roof function that is a weak solution satisfying the boundary conditions 
almost everywhere. Thus, this provides a pathological example of an exceptional 
domain in a weak sense. In order to construct such an unbounded non-Smirnov 
domain, let us follow the method in the above mentioned [13 , which is presented 
in Duren's book [T2J Section 10.4]. We recall that the construction is carried out 
by "working backwards", first writing down a singular inner function S(z) as a 
candidate for the derivative f'(z) of the conformal map f(z). The difficulty is then 
to show that f(z) is not only analytic, but is also univalent so that it actually gives 
a conformal map from D to some domain O. Univalence is established using a crite- 
rion of Nehari which states that the following growth condition on the Schwarzian 
derivative (Sf)(z) is sufficient for univalence: 

(2.2) (Sf)(z) < 

Let us follow this procedure, indicating the step that needs to be modified. Start 
with a measure [i < 0, singular with respect to Lebesgue measure on the circle, yet 
sufficiently smooth, so that it belongs to the Zygmund class A* (cf. [HI Section 
10.4]). 

We will also require [i to have the first moment zero, i.e., 

t-2-K 

(2.3) / e ie d^{9) = 0. 

Jo 

This can always be acheived by symmetrizing fj, around the origin and re placi ng fi 
by |(<i/x(0) + dfi(—6)). Then the center of mass is at the origin, which is (2.3). 
As in [12 , let F(z) be the Schwarz integral of fj, 

1 C 2n e ie + z 

Let g(z) be the exponential of a constant (to be chosen later) times F, 

g(z) = exp{-aF(z)}. 

Here is where we depart slightly from |12j in order to get an unbounded domain as 
the image of f(z). Instead of taking g(z) as a candidate for f'(z), we take 

f{z)=g(z)/z 2 . 

Note that the residue of f'(z) is zero (from having made the first moment of [i zero) 
so that its antiderivative f(z) is analytic in D except for a simple pole at z = 0. 
Also, \ f'(z)\ = 1 a.e. on dD. 

A calculation shows that the Schwarzian derivative Sf of / is: 

(Sf)(z) = (Sg)(z) - = -aF"(z) ^F'{zf + ^F'(z). 

z g\zj l z 

As explicitly stated in [TJ Section 10.4], F"(z), F'{z) 2 , and F'(z) are each 

°((T^N) 2 , 
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Moreover, by the vanishing of the first moment of fi, F'(0) — 0, so that F'(z)/z is 
also O ( (i_}z|)2 ) • Thus, for a small enough choice of a, (Sf)(z) satisfies the Nehari 



criterion for uni valence (2.2). 



Hence, f(z) is a conformal map mapping {|z| < 1} onto the complement of a 
Jordan domain with rectifiable boundary. To see why the boundary is rectifiabile, 
note that, as stated in [321 Section 10.4]), g(z) e H , and so f'(z) — g(z)/z 2 is in 
H 1 in an annulus < r < \z\ < 1. 

This seemingly excessive construction of an exterior pseudocircle cannot be 
avoided by simply taking an inversion of an interior pseudocircle; the result will 
be non-Smirnov, but it will not be an exterior pseudocircle. Nor can one simply 
take the complement. As P. Jones and S. Smirnov proved in [22], the complement 
of a non-Smirnov domain is often Smirnov! (This unexpected resolution of a long 
standing problem put to rest all hopes to characterize the Smirnov property in 
terms of a boundary curve.) 

Remark 2. The above examples of non-Smirnov exceptional domains lead to as- 
suming fl is Smirnov in our main theorems (but we allow u to be a weak solution). 

An alternative approach is to require u to be a "classical solution" that satisfies 
the boundary condition everywhere (and not just almost everywhere), then non- 
Smirnov domains are ruled out. Moreover, real-analyticity of the boundary then 
follows automatically. To be precise, we have the following Lemma. 

Lemma 2.2. Iffl C M 2 is exceptional and the roof function u is a "classical solution" 
in C 1 (fl), then dfl is locally real- analytic. 

Proof. The analytic completion f(z) = u + iv (possibly multivalued) maps fl into 
the right halfplane, since u is positive. The Neumann condition for u implies that 
\f'(z)\ = 1 on dfl. Also, u G C 1 ^) implies that /' G C(H). 

Choose a point zq G dfl, and let = f(zo). Let g(Q = / _1 (C) denote the 
local inverse of f(z). Choose a neighborhood U of £ and let F := U fl {5R(C) > 0}. 
Choose U small enough so that g G C(F). 

Since |g'(C)| = 1 on dfl, we can also choose U small enough that g' does not 
vanish in F. This implies that h(Q = Log(</(£)) is analytic in the interior of 
F and continuous in F. We have 5i{ft.(£)} vanishes on the imaginary axis, since 
|ff'(C)| = 1 there. Thus h(Q extends to a neighborhood of £o by the Schwarz 
reflection principle. This allows us to extend g'{z) and therefore g{z) and f(z) 
extend analytically across Zo, since u := 5R/ = on dQ and |Vtt| = 1 on dfl near 
zq. The lemma is proved. □ 

Corollary 2.3. If dfl is C 2 -smooth and fl is exceptional then dfl is locally real- 
analytic. 

Proof. C 2 -smoothness of dfl implies that u is in C 1 (f2). It remains now to refer to 
Lemma 2.2. □ 

Using Kellogg's theorem on regularity of conformal maps up to the boundary, cf. 
[3T1 Ch. 3], one easily extends the above corollary to C 1,Q , a > 0, boundaries and 
even merely to C 1 boundaries. We shall not pursue these details here. It would 
be interesting to find sharp necessary and sufficient conditions for the a priori 
regularity of the boundary that would guarantee the conclusion of Corollary |2.3| 
As we have mentioned in the beginning of this section, it is necessary to assume that 
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the domain is Smirnov, but it is not at all obvious that this is indeed sufficient, cf. 
a related discussion in [TU] and [TT] regarding nonconstant functions in E p classes 
with real boundary values. 

3. The Case When Infinity is an Isolated Boundary Point 

Theorem 3.1. Suppose CI is an exceptional domain whose complement C \ CI is 
bounded and connected, and assume CI is Smirnov. Then CI is the exterior of a 
disk. 

Proof. Let u be a roof function for CI. Positivity of u implies, by Bocher's Theorem 
[21 Ch. 3], u(z) = Uq(z) + Clog \z\ for some constant C, where uq{z) is harmonic 
in CI U {oo}, and uq(z) approaches a constant at infinity (the "Robin constant" of 
dCl). Thus, in view of the Dirichlet data of u, u(z) is a multiple of the Green's 
function of CI with a pole at infinity, and taking v(z) to be the harmonic conjugate 
of u{z)/C, we have a conformal map g(z) — e u ( z )/ c + lv ( z ) from SI to the exterior of 
the unit disk (note that g{z) is single- valued in SI). 

Using both the Dirichlet and Neumann data, we have |</(z)| = 1/C a.c. on <9S1, 
and therefore 

a.e. on <9B. Since g 1 has a simple pole at infinity, (g 1 )' is analytic. Also, (g 1 )' 
is in iJ 1 (C \ D) since dCl is rectifiable. Since SI is Smirnov, the latter function is 
outer and also has constant modulus on the unit circle a.e., which together imply 
that it is constant. (Recall from Section [2] that by formula (2.1 1 an outer function 



is determined from its boundary values.) Hence g 1 is a linear function and dCl is 
a circle. □ 

We defer proving a higher-dimensional version of this result until Section [7J but 
we mention here that under more smoothness assumptions the higher-dimensional 
case can be proved using a theorem of W. Reichel [32]. 

Under additional smoothness assumptions, the hypothesis of Theorem [3T] guar- 
antees that CI is a special type of arclength quadrature domain. The following is 
then an immediate corollary of a result of B. Gustafsson [TH1 Remark 6.1]. 

Theorem 3.2 (B. Gustafsson, 1987). Suppose CI is a finite genus exceptional domain, 
with piecewise-C 1 boundary, and infinity is not a point on the boundary of CI. Then 
CI is the exterior of a disk. 

This removes the condition that the complement of SI is connected. 

Proof. We will show that SI is an arclength null quadrature domain for analytic 
functions vanishing at infinity. At first, consider as a class of test functions to 
integrate over dCl rational functions r(z) in CI vanishing at infinity. 

Let f{z) = u(z) + iv(z) be the analytic completion of the roof function u. Note 
that f'(z) is single- valued (since it is the conjugate of the gradient of u), and by 
Bocher's theorem cited above f'(z) = 0(|z| _1 ). Since the gradient of u is the 
inward normal of <9S1, f'{z) — jtj^ is as well. The unit tangent vector ^ is a 
90-degree rotation of the normal vector jT^y- Thus, if'(z)dz = ds. We then have 
a quadrature formula for integration of r{z) with respect to arclength: 

(3.1) / r(z)ds = i I r(z)f'{z)dz = 0, 

Jan Jan 
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where the vanishing of this integral is obtained by deforming the contour to infinity 
where f'(z)r(z) — 0(|z|~ 2 ). Indeed, r(z) = 0(|z| _1 ) by assumption on the test 
class, and f'(z) = 0(|z| _1 ) as mentioned above. 

If the boundary of is piecewise-C 1 then the rational functions are dense in E p 
classes (see [121 Thm. 10.7], and for the multiply connected case, [39], [40], [41]). In 
particular, the rational functions r(z), vanishing at infinity, are dense in the space 



of functions E(Q,) considered in [IB]. Thus, (3.1| shows that fl is an arclength 
null quadrature domain for this space of functions, and the result now follows from 
Remark 6.1 in [TS]. □ 

4. The Schwarz Function of an Exceptional Domain 

The Schwarz function of a real-analytic curve T is the (unique and guaranteed to 
exist near T) complex-analytic function that coincides with z on T. For the basics 
on the Schwarz function we refer to [9] and [56] , 

We recall two basic facts needed in the proof of the next proposition. 

(i) OnT, \S'(z)\ = l. 

(ii) The complex conjugate of \f—S'(z) is normal to T. 

Statement (i) follows from the chain rule and the fact that the complex conjugate 
of the Schwarz function, S(z), is an involution (see [5] Ch. 7]). Statement (ii) 
follows from the formula for the complex unit tangent vector T(z) = = , 1 

1 6 W da y/Wlz) 

expressing the derivative of z with respect to the arc-length along T (see again 
Ch. 7, Formula (7.5)]). 

Proposition 4.1. If £1 is an exceptional domain such that the roof function is a 
classical solution, then the z-derivative of the roof function is given by u z (z) — 
CyJ—S'{z), where c is a real constant and S(z) is the Schwarz function of d£l. In 
particular, S' (z) is analytic throughout £1. 

Remark 3. If, for instance, the constant Neumann data for the roof function is 1, 
then the constant above c = ±1/2 where the sign depends on the orientation of the 
boundary. 



Proof. Lemma 2.2 implies that L is locally real-analytic. So T has a Schwarz func- 
tion S(z). The complex conjugate of the analytic function u z is normal to T (since 
u has zero Dirichlet data). In light of the constant Neumann data, we then have 



|u z (z)| = |(m 2 (z))*| = \\u x + iu y \ — \Jv%. + u 2 is constant on T. This, along with 

the statements (i) and (ii) above, shows that on T the vectors u z (z) and yf-S'(z) 
are parallel and each have constant length. Therefore, for some real constant c, 
the equation u z (z) = cy/—S'{z) holds on V. But since u z and y/-S'(z) are both 
analytic, the equation is true everywhere that either side is defined. In particular, 
this guarantees analytic continuation of S'(z) throughout il. □ 

Let us use the Schwarz function to give a heuristic argument that the boundary 
of an exceptional domain can pass through infinity at most twice. In fact, the angle 
between consecutive arcs at infinity must be ir (and obviously there cannot be more 
than two such angles at infinity). Suppose the boundary of a domain has a corner 
where two arcs meet at an angle different from 0, tt, or 2ir. Then the derivatives 
of the Schwarz functions of the two arcs have a branch cut along a third arc that 
propagates into the domain from the corner. To see why this is the case, note that 
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the Schwarz function of an arc can be approximated near a point by the Schwarz 
function of the tangent line. Thus, to first order, the jump along the branch cut 
is linear, so to zeroth order, the jump of S' is determined by the slopes of the 
tangent lines. If the angle is or 2ir then the tangent line is the same for each 
arc, but the orientation changes, so there is still a jump due to the sign change. In 
the case of an angle of ir both the tangent line and the orientation are unchanged. 
Thus, for any angle other than 7r, S'(z) has a jump across a branch cut between 
the two boundary boundary components. For an exceptional domain, u is a global 
solution throughout fl, and so Proposition |4. 1| indicates that the Schwarz function 
cannot have such branch cuts. Thus, the angle between consecutive boundary arcs 
at infinity can only be 7r, and there can be at most two such angles. 

In the above informal argument, we have assumed that each arc is real analytic 
at infinity, so that the Schwarz function has an expansion. Alexandre Eremenko 
related to us the following proof [TS] using ideas from [5] that extend techniques 
due to Ch. Pommerenke. No regularity assumptions on dfl are required. Also, an 
important part of the theorem readily extends to higher dimensions. 

We recall that a Martin function is a positive harmonic function M in a region 
n with the property that for any positive harmonic function v in f2, the condition 
v < M implies that v — cM, where c > is a constant. (Often, Martin functions 
are called minimal harmonic functions - cf. [20J .) Martin functions on finitely 
connected domains are simply Poisson kernels evaluated at points of the Martin 
boundary, the boundary under Caratheodory compactification (prime ends) of the 
domain (see [5]). 

Theorem 4.2 (A. Eremenko [15 ). The roof function u of any exceptional domain £1 
is a convex combination of at most two Martin functions of at infinity. Moreover, 
u(z) — 0(\z\), and in two dimensions we also have Vu(^) = 0(1) in fl. 

Remark 4. M. Traizet obtained the estimate |Vu| < 1 in fl for domains with 
finitely many boundary components using the Phragmen-Lindelof principle. For 
Smirnov domains fl it suffices to show that u z belongs to the class N + (cf. [10] ) 
in order to conclude that the analytic function u z is bounded by 1 in fl. However, 
even this assumption is not needed here, and it is possible to establish the estimate 
on Vu in full generality. Alexandre Eremenko has kindly permitted us to include 
his argument here. 

Proof. First we note that, as observed in [51 Lemma 1], if u is a positive harmonic 
function in a disk (or a ball in higher dimensions), D(a, R), of radius R centered at 
a, and u{z\) = for some boundary point z%, then 

(4.1) u(a) < 2R\Vu( Zl )\. 

This immediately follows from Harnack's inequality for D(a,R) as for z E D(a,R) 
u(a) u{z) u(z) — u{zi) 



R+\z — a\ 1 — \z — a\ 1 — \z — a\ 



and letting z — > Z\ establishes (4.1 1 



Applying ( |4.1[ ) when a € CI and R is the distance from a to dfl, gives u(a) < 
2R < 2(\a\ + const). So u{z) = 0(\z\), as z -> oo. 
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The fact that u is a combination of at most two Martin functions now follows 
from a standard argument using Carleman's inequality, see for example [57]. For 
the higher-dimensional case, one must use [TTJ instead of |27j . 

Next we show, in the two dimensional case, the additional claim that Vm(z) = 
0(1). Let R > and consider an auxilliary function 

|V«| 
u + R 

where R > is a parameter. A direct computation shows that 

(4.2) Alogu^>«4 
and wr(z) = 1/R for z € dD. We claim that 

(4.3) w R (z) < 2/R, zeD, 

from which the result follows by letting R — > oo which gives |Vu| < 2 in f2. 



Suppose, contrary to (4.3), that w_r(zq) > 2/R, for some zq G ^- Let 



2_R 

u ( z ) = ^2 i 12 ' z e ^(^O' -R) = : l z - z o\ < R}- 

R 1 - \z — zqY 

Obviously, v(z) > 2/R. A computation reveals that Alogt* = v 2 . Let 

K = {z G H n D(z 0) -R) : w fl (z) > w(z)}. 

We have Zq € if, since t>(zo) = 2/i2. Let K be the component of K, containing z - 
Then we have wr{z) = v(z) on 8Kq, since wr(z) < v(z) on 9f2 n D(zq,R) while 
v(z) = +00 on 8D(zq, R). On the other hand, 

A(\og w R - log v) > wjj - v 2 > in ivT - 

So the subharmonic function log u — log v is positive in Kq and vanishes on the 
boundary, a contradiction. □ 

Remark 5. This a priori estimate implies the following corollary showing that the 
boundaries of exceptional domains are extremely regular. Namely, they are locally 
real analytic and even parameterized from the unit circle by the antiderivative of 
a rational function. In particular, it validates the preceding argument using the 
Schwarz function, and establishes that the boundary passes at most twice through 
infinity each time with an angle of tt. The only additional assumptions needed here 
are that the domain is Smirnov (cf. Section 2) and simply connected. 

Corollary 4.3. Let Q be a simply connected Smirnov domain, and let h(C) be the 
conformal map from D to fi. If Q is exceptional then is a rational function, 

and either: 

Case (i). hi has one pole on dD, or 

Case (ii). h! has two poles on 3D. 

Proof. Let u be a roof function for f2, and f(z) = u + iv its analytic completion. 
Since u > 0, f(z) takes f2 into the right halfplane, and f(h(()) takes the unit disk 
D into the right halfplane. Adding an imaginary constant if necessary, we may 
assume that f(h(0)) > is real. Then, by the Herglotz Theorem (see [2TJ Ch. 3], 
[L2l Ch. 1]), we can represent /(/i(C)) as 

(4-4) /(MO) - / e ir^d»(9), 



OVERDETERMINED PROBLEM IN POTENTIAL THEORY 



11 



with /i positive. 

Now since is the pull back to D of the roof function u, which by Theorem 

|4.2| is a convex combination of at most two Martin functions, [i consists of at most 
two atoms. 



Thus, differentiating (4.4) 



(4.5) f(h(0)-h'(0=R(0, 

where R(() is a rational function with either one or two double poles on dD (at the 
atoms of (tx). Since f'(h(()) is a bounded analytic function in D with \f'(h(Q)\ = 1 
a.e. on dil, f'(h(()) is an inner function. Moreover, h'(Q is an outer function, 
since D, is Smirnov. 

For a rational function such as R(C) the canonical factorization given by Theorem 
12. II reduces to: 

R(C) = B(0 ■ F(C), 

with B a Blaschke product and F a (rational) outer function. (The singular factor 
S(() is trivial, since R(() has no essential singularities.) By the uniqueness of the 
canonical factorization, h'(() and f'(h(£)) equal F(Q and -B(C) respectively (up 
to multiplication by a unimodular constant). Hence, h'(() = F(() is rational, and 
f(h(C)) = B(z) is a Blaschke product. □ 

5. The Case When Infinity is a Single Point on the Boundary 

In this Section, we characterize the halfplane as the only simply connected ex- 
ceptional domain having infinity as a single point on the boundary. This extends 
[19l Prop. 6.1] by removing the additional hypothesis d x u > in [19] on the roof 
function in fi. 



Theorem 5.1. Suppose Q satisfies the hypothesis of Corollary 4-3 Then Case (i) 
implies that Q is a halfplane. 

Remark 6. (i) We have not been able to drop the assumption that £1 is simply 
connected, but as mentioned in the introduction and Section 9, M. Traizet has 
recently established the result under the assumption of finitely many boundary 
components. |38j . We note that, in the simply connected case, this assumption is 
stronger than ours (since we allowed for infinitely many boundary components in 
Section 4). 



(ii) We have not been able to prove a higher dimensional version of Theorem 5.1 
(cf. Section 7). 

Proof. Using the same notation / and h from the proof of Corollary |4.3[ 



/(MC))= / 3^4*0 

Job e ~ t> 

for some finite positive measure ji on dH>. By assumption, we are in the case when 



h! has one pole, and according to the proof of Corollary 4.3 /i is an atomic measure 
with a single point mass. Without loss of generality, we can place it at the point 
e lf> = 1. 
Thus, 

(5T) /(ft ( C )) = ci±|. 
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Differentiating (5.1), 
(5.2) 



f'(h(0)h'(C) = 2C 



(i-C) 5 



and as asserted in the proof of Corollary 4.3 f'(h(C)) is the Blaschke factor of the 
right hand side, which has no zeros, so f'{h(Q) is a unimodular constant. Therefore, 
/ is a linear function and is a halfplane. □ 



6. The Case When Infinity is a Double Point of the Boundary 

In this section we characterize the nontrivial example found in [19 . Suppose f2 is 
a simply connected domain and f2 is exceptional. By Corollary |4.3[ recall that the 
derivative ft-'(C) of the conformal map from the disk onto f2 is a rational function 
with either one or two double poles on <9B. 



Theorem 6.1. Suppose fi satisfies the hypothesis of Corollary 4-3 Then Case (ii) 
implies that is, up to similarity, the image of the strip |3u>| < it/ 2 under the 
conformal map g(w) = W + sinh(w), while the analytic completion of the function 
u(g(w)) is the function f(g(w)) — cosh(w). 

Remark 7. The exceptional domain Q that is the image of the strip under the con- 
formal map C — > C + sinh £ is precisely the exceptional domain found by the authors 



Theorem 6.1 together with Theorems 3.1 and 5.1 and Theorem 4.2 show 



that this f2 is, under an assumption on the topology essentially the only nontrivial 
example of an exceptional domain in K 2 . It turns out that some assumption on 
topology is necessary as there is yet a whole one-parameter family of non-similar 
exceptional domains that have infinite genus (See Section 9). However, under the 
assumption of finitely many boundary components, the example described in The- 



6.1 is the only nontrivial example as the previously mentioned recent work of 



M. Traizet shows 138 . 



4.3 



and Theorem 



Proof. Using the same notation as in the proofs of Corollary 
5.1 we have that h'(Q is a rational function, and according to (4.5| f'(h(()) is as 
well. This justifies applying the argument principle to study f(h(Q) and /'(/i(C))- 
Namely, we will prove the following. 

Claim: The function / solves a differential equation: 



(6.1) 




z e fi, 



after simple normalizations described below. 

Before proving this Claim we solve the differential equation to see that it gives 
the desired result. Separating variables, 



l f + l 
/-I 



df = z + C. 



Make the substitution / = cosh(u>),z = w 
gration C = 0). Now using the conditions 



sinh(w) (fixing the constant of inte- 



ft/(z(iu)) = for z £ dQ, and $tf(z(w)) > for ( e fl, 



OVERDETERMINED PROBLEM IN POTENTIAL THEORY 



13 



and the identity 9?cosh(a: + iy) = cosh(x) cos(y), we find that the pre-image of the 
domain in the w-plane is the strip \Qw\ < n/2. Therefore, SI can be described as 
the image of the strip under the map z(w) — w + sinh(w). 

□ 

Proof of Claim. We will use the argument principle to show that both sides of the 



Equation (6.1 1 provide a conformal map from SI to D. 

Starting from the formula which relates the tangent vector T(z) on <9S1 and the 
derivative of the analytic completion / of u(z), 

(h = -% = 1 

(z> ds f'(z) ^wjry 

we obtain from the continuity of T(z) through the double point at infinity (see 
Figure 1), that 

I d\ogf{z) = 2m. 
Jan 

We conclude that /' is a single-sheeted covering of the unit disk by the domain SI, 
and that it has only one zero, at some point zq 6 Q. 

We may assume that f(zo) = 1. If not, say f(z ) = a + ib 7 a > 0, then one 
may subtract the constant ib from / (this just amounts to choosing a different 
harmonic conjugate for the same roof function), so we have f(z ) = a. Then one 
may simply replace 1 with a in the claim, and integrating the differential equation 
is done similarly resulting in a dilation of the original solution. 

Consider now the function defined on SI, taking values in the unit disk D 



' /(*)-! 



Then g is also a univalent map from SI into D. Indeed, by the argument principle, 
jfefej; is a branched, two-sheeted covering of the disk, since it maps each of the 
two boundary components shown in Figure 1 onto T, Moreover, the single branch 
point zq is mapped to the origin, so that taking the square root gives a single- valued 
analytic function. 

Also, f'(zo) — g{zo) — 0. This uniquely determines the conformal map up to 
a unimodular constant, which we may assume is 1 (after a rotation), and we then 



arrive at the differential equation (6.1) 



7. An Extension of Theorem 13.11 to Higher Dimensions 



□ 



In this section, we notice that some results in Section [3] extend to higher dimen- 
sions. 

Theorem 7.1. Suppose SI is an exceptional domain in W whose exterior is bounded 
and connected. If 90 is C 2 ' a -smooth, a > 0, then dfl is a sphere. 

Proof. Let u be a roof function for SI, and let v(s) = | s ^_ 2 denote the Newtonian 
kernel. Fix y £ fl and take a small ball B e centered at y. Take also a large ball Br 
of radius R that containes both B e and the complement of SI. 
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Figure 1 . Local geometry of the boundary T = dil near infinity. 



Since u(x) and v(x — y) are harmonic in Q, \ B e , Green's second identity gives 



(7.1) / (v(x - y)d n u(x) - u(x)d n v{x - y)) da x = 0. 

JdB R +dn-dB s 

Letting R —¥ oo, we can drop the integration over OBr, since again by Bocher's 
Theorem [31 Ch. 3], near infinity u(x) s» |a;| 2 ~™. 
Since, u{x) = on dfl and d n u(x) = 1 on 9f2, 



(7.2) / v(x - y)dcj x = I (v(x — y)d n u(x) — u(x)d n v(x - y)) da x . 

Jon JdB e 

Let U be the bounded domain such that W l \ U = Q. The outward normal for 
dU is opposite to that of dil, and since v(x — y) — -^i on dB e , 



(7.3) / v(x-y)da x =l {-j^2d n u(x) + u(x)d n v(x-y))da x . 

JdU 



For the first term on the right-hand-side, we have 

/ 7^=2 d n ii(x) da x = / Au(x)dV = 0. 

JOB* JB r 



So, 



v(x - y)da x = / u(x)d n v{x - y)da x -> u(y), 

dU JdB E 



as e — > 0. So, u(y) is the single layer potential with charge density one on the surface 
dU. That U is a ball now follows from a theorem of W. Reichel [3H Theorem 1]. 

□ 

Remark 8. Reichel's result holds for more general elliptic operators than the Lapla- 
cian. In the setting of the Laplacian, J. L. Lewis and A. Vogel [55] characterized 
the sphere in terms of its interior Greeen's function under weaker regularity as- 
sumptions, namely, the boundary is assumed Lipschitz. In that case, the Neumann 
condition can be assumed to hold almost everywhere on the boundary. Thus, the 
hypothesis of Theorem 1 7 . 1 1 could be weakened by checking that the same proof [23J 
works for the exterior case we are interested in. Yet, we have chosen an easier and 



OVERDETERMINED PROBLEM IN POTENTIAL THEORY 



15 



more transpoarent path to apply Rcichel's result directly, even though it requires 
a stronger regularity on the boundary. 



8. Nonexistence of a Higher-dimensional Analog of the cosh(z) 

Example 

The authors in [115] expressed a suspicion (see Remark 2.1 in [TIJJ) that there ex- 
ist n-dimensional, rotationally-symmetric examples similar to the two-dimensional 
example {(x, y) E R 2 : \y\ < f + cosh(x)} that appeared in Section [6[ We show 
that there does not exist an exceptional domain in R whose boundary is generated 
by rotation about the cc-axis of the (two-dimensional) graph of an even function. 

Theorem 8.1. There does not exist a rotationally-symmetric exceptional domain O 
in R 4 that contains its own axis of symmetry and whose boundary is obtained by 
rotating the (two-dimensional) graph of an even real-analytic function about the 



Remark 9. (i) Our proof will rely heavily on two tricks, one exploiting the assump- 
tion that n = 4, and the other using the assumption that the generating curve 
is symmetric. However, we strongly suspect a more general non-existence of such 
examples in R™ for any n > 2. Therefore, we conjecture the following. 

Conjecture 8.2. For n > 2, there does not exist an axially symmetric, exceptional 
domain in R™ that contains its own axis of symmetry. 

(ii) The assumption that the domain contains its axis of symmetry rules out 
the exteriors of balls and circular (or spherical) cylinders, respectively (which are 
clearly exceptional domains as was noted in [EI]). Also, A. Petrosyan and K. 
Ramachandran pointed out to us that the nonconvex component of the exterior of 
a certain cone is also an exceptional domain. In R 4 , using the a:-axis as the axis of 
rotation, the cone is the rotation of {(x,y) : y 2 — x 2 = 0}, and the roof function 
in the meridian coordinates x,y where y is the distance to the a:-axis in R 4 , is 
u{x,y) = v ~ x for y > 0. 



Proof of Theorem 8.1 Suppose that il is such a domain in R . Namely, the bound- 
ary dQ, is obtained from rotation of 7 := {(x,y) € R 2 : y = g{x)}, with g{— x) = 
g{x). i.e., the boundary of O is given by 

{(xi,x 2 ,x 3 ,x 4: ) e R 4 : \Jx\ + x\ + x\ = g{xi)}. 

Considering the boundary data, the rotational symmetry of the domain will be 
passed to the roof function, so that, abusing notation, we can write 

u(xi,x 2 ,x 3 ,x 4 ) = u(x,y). 

For clarity, we emphasize that the x-axis corresponds to the axis of symmetry and 
the y-coordinate gives the distance from the axis of symmetry. 

For axially symmetric potentials v in R™ the cylindrical reduction of Laplace's 
equation is: 

. (n — 2)v v 
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= 0, if and only if yu(x,y) is a 
yAu + 2u„ 



where x = X\ and y = \fx\ + ... + x\. Moreover, in the case we are considering, 
when n = 4, u satisfies the equation Au + 
harmonic function of two variables x and y. Indeed 

A(yu) = yAu + 2Vu • Vy + uAy — y^u, -r ^u, y 

(The trick that reduces axially symmetric potentials in ]R 4 to harmonic functions 
in the meridian plane is well known, cf. [26 and [23 .) 

Since yu(x, y) is then harmonic in the unbounded two-dimensional domain D 
bounded by 7 and its reflection (which we denote by 7) with respect to the x-axis, 
this implies (yu(x,y)) is analytic in the domain D, where as usual z — x + iy. 
The Cauchy data (originally posed in R 4 ) imply that u z = \(u x — iu y ) coincides 
with y/-S'(z) on 7 and 7. This implies that the analytic function 



(8.1) 



W{z) 



(y u ) z = y u 



yu z 



coincides with z yJ—S'(z) on 7 and 7. Since z sj— S'(z) is analytic, this 
actually gives a formula for W(z) valid not only on 7 and 7: 



(8.2) 



2/ 



We note that (8.2) can be used to analytically continue S(z) to all of D, but this 
is not needed in our proof. 

Let /(C) be the conformal map from the strip S := < ^} to Z? such that 

/(0) = and arg{/'(0)} = 0. The two-fold symmetry of D implies that /(C) is an 
odd function. Indeed, otherwise h(£) = — /(— C) gives another conformal map from 
the strip E to D. But, h(0) = -/(0) = and ft'(Q) = /'(0) implies h = f, by the 
uniqueness of the conformal map (up to choice of /(0) and argument of /'(0)). 

The Schwarz functions St, Sb of the top and bottom edges of the strip E are 
"St(0 = C~ h an d 5*b(C) — C + *- I n terms of the conformal map /(£), the pull-back 
to the C-plane of the Schwarz functions S+ and S- of 7 and 7 (resp.) satisfy (see 
Ch. 8, Eq. 8.7]) 



(8.3) 
(8.4) 



S±(/(0) = /(C=Fi), and 

/'(CtO 



S±(/(0) 



/'(C) 



Substituting these into (8.2 1, we obtain two expressions for the pullback of W(z) 
to the strip E: 



(8.5) 



/(C) -/(CtQ / / ; (CTi) 
2* V /'(C) 

Even though M^(/(C)) is analytic throughout E, we caution that these two ex- 
pressions (one expression for "+" and one for "— ") may only be valid near the 
bottom and top sides (respectively) of the strip E. 
Claim: The function W(f(Q) is odd. 

Before proving the Claim, let us see how it is used to finish the proof of the 



Theorem. The fact that W{f(()) is odd implies W(0) = W(f(0)) = 0. By (8.1 1 we 
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then have that -J-u + yu z vanishes at z — 0, which implies that w(0,0) = 0. This 
contradicts the positivity of u. 

Proof of Claim. We wish to show that V(C) = W(f(Q) + W(f (-(,)) vanishes 
identically. Use each of the expressions in ( |8.5[ ) above to represent W(f(£)) and 
W(f(— ()), respectively. 



r»tf vin /(0-/(C-») f(C-0 , f(-()-f(-( + i) fK+j 

(8 ' 6) no= 2i V^^cT a r~JFcr 

We show that this formula vanishes where it is valid, which then implies that 
V(C) vanishes identically throughout S. For this, we use the fact that / is odd and 
consequently /' is even. 



(8J) nc) - * V~~ra~ k V~~7(0~ 

This establishes the Claim. □ 

□ 

9. Concluding Remarks and Main Conjecture 

1. It is tempting to conjecture that the three examples in the plane studied 
above are the only exceptional domains in the plane as suggested in [19] . 

However, there is a remarkable family of infinitely-connected exceptional do- 
mains. These were discovered as solutions to a fluid dynamics problem by Baker, 
Saffman, and Sheffield in 1976 0]. See also [5] for a more detailed account. The 
original problem there was to find hollow vortex equilibria with an infinite periodic 
array of vortices, i.e., "spinning bubbles" amid a stationary flow of ideal fluid. The 
domain occupied by fluid turns out to be an exceptional domain with an infinite 
periodic array of holes, and the roof function is a stream function of the fluid flow, 
see Figure [3j The constant Dirichlet condition corresponds to the requirement that 
the boundary of each hollow vortex is a stream line, and the constant Neumann 
condition corresponds to the requirement that the fluid pressure should be balanced 
at the interface by the pressure inside each bubble which is assumed constant. The 
latter correspondence is more subtle; in order to have constant pressure along a 
stream line, the fluid velocity (which equals the normal derivative of stream func- 
tion) should be constant according to Bernoulli's law. 

This infinite genus example leads us to add to the conjecture the assumption 
that the domain has finite genus. 

Conjecture 9.1. The only finite genus exceptional domains in M 2 are the exterior 
of the unit disk, the halfplane, and the domain described in Theorem \6.1\ 



Remark 10. As mentioned in the introduction, Martin Traizet |38j recently an- 
nounced a classification of exceptional domains. His results confirm our conjecture 
for domains having finitely many boundary components and also show that the 
above infinite genus example is the only periodic exceptional domain for which the 
quotient by the period has finitely many boundary components. His methods use 
a remarkable nontrivial correspondence to minimal surfaces, perturbing an excep- 
tional domain by harmonically mapping it to another domain in such a way that 
the graph of the new height function (which pulls back to the roof function in the 



18 



DMITRY KHAVINSON, ERIK LUNDBERG, AND RAZVAN TEODORESCU 




Figure 2. An infinite- genus exceptional domain that also provides 
a hollow vortex equilibrium. Level curves of the roof function are 
stream lines. The shape of the bubbles ensures that the pressure 
dictated by Bernoulli's law is constant at the fluid-bubble interface. 




Figure 3. There is actually a whole one parameter family of dif- 
ferent bubble shapes. As noticed in 38 , each of the three previ- 
ously known examples can be recovered as different scaling limits 
of this family. In that sense, this family includes all known exam- 
ples. 

original domain) satisfies the minimal surface equation. A miraculous (and crucial 
to his approach) by-product is that, whereas the graph of the roof function meets 
its boundary at a 45-degree angle, the minimal graph meets its boundary vertically 
so that gluing it to its own reflection over the sy-plane results in a smooth minimal 
surface (without boundary!) embedded in E 3 . 

2. Regarding the higher-dimensional case, we conjecture the following extension 
of Theorem |5.1 1 to higher dimensions. 

Conjecture 9.2. Suppose f2 is an exceptional domain in R" that is homeomorphic 
to a halfspace. Then 17 is a halfspace. 

3. The connection to the Schwarz function in Section [4] reveals that exceptional 
domains are arclength null- quadrature domains. That is, for any function /, say 
analytic in fi, continous in VL, integrable over the boundary, and decaying suffi- 
ciently at infinity, we have fgafds = 0. Indeed, Jg^fds = j dn f(z)j^T(z)ds — 

Ion f( z )\/S'( z )dz, where T(z) is the complex unit tangent vector (see Section jSj), 
and now this integral vanishes as long as the integrand decays sufficiently at infinity. 
Null-quadrature were previously studied in the case of area measure. They were 



OVERDETERMINED PROBLEM IN POTENTIAL THEORY 



19 



characterized in the plane by M. Sakai [33] - Our current study can be seen as a 
step toward characterizing null-quadrature domains for arclength. 

4. Other interesting connections involve differentials on Riemann surfaces. The 
study of Gustafsson [T5] used half-order differentials on the Schottky double of an 
arclength quadrature domain. From a different point of view, the boundary of 
an exceptional domain is a trajectory of the positive quadratic differential — (<i/) 2 , 
where f(z) is the analytic completion of the roof function. 



5. The differential equation (6.1| can be solved by a more general substitution us- 
ing Jacobi elliptic functions [TJ p. 567, §16]: /(£, k) = cos(#)cn(£, k) + sin(6')sn(C, k) 

and z(C) = 4>(C) + cos(6»)sn(C, k) - sin((9)cn(C, k), where 0(C) = f dn(f , k)d£ and 6 
is an arbitrary phase, 6 € [0, 2ir}. 

For a given value of the elliptic modulus k € [0, 1], we define the corresponding 
domain F through its fundamental periods T\{k) — AF(-k/2, y/1 — k 2 ) and T^/c) = 
4F(7r/2, k), where F(w/2,k) = K(k) is the complete elliptic integral of the first 
kind [1] p. 590, §17.3]: 

«r/2 -, 

K(k)= / =d6 
Jo y4 - k 2 sm 2 (6) 

It diverges for k = 1 and equals ir/2 for k — 0. 



Then it is straightforward to check that (6.1) is satisfied by f(z), due to the 
identity [TJ p. 573, §16.9]: 

1 = [-sn(z) cos(0) + cn(z) sin(6»)] 2 + [cn(z) cos(6») + sn(z) sm(9)} 2 . 

Let 7 be the pre-image of dfl under z(£): it consists of two pieces 7±>7- = —7+, 
dividing the fundamental domain F into three sub-domains. Denote the component 
which contains the origin by D , then since /(0) = 1, we conclude that $if(z) > 
for z € Dq \ 7±, and we have proven the following result. 

Proposition 9.3. The exceptional domain f2 is the image of the domain Do(k) under 
the map z(Q = <fi(() + cos(9)sn(C t , k) — sin(0) cn(C, , k). 

Remark 11. The case discussed in the proof of the theorem corresponds to the 
degenerate elliptic modulus k = 0. Then the domain F becomes the infinite strip 

T'i(O) = 4:K(1) — y 00, T 2 (0) = 4K(0) = 2tt, 

while the functions /, g become (using the fact that dn(z, 1) = 1) 

z(C) = C + sinh(C), f(z(0) = cosh(C). 

As noted before, the conditions !ft/(()| 7± = give the pre-image j± := z^ 1 (dfl) — 
{3£ = ±|}, and the pre-image of the domain, D , becomes the strip \Ss(\ < J. 

6. Note that the domain D$(k) is the pre-image of the unit disk under the map 
C(io) : F -> D, 

> , . sa(w, k) — i , .„ . 

CM= ' , 0,1, 

sn(w, k) + 1 

with the support of /i at points C± = iyq^ff , where /i is the measure discussed in 
the proof of Corollary |4.3[ The case k — > corresponds to the strip domain and 
to (± = ±1. The reparametrization invariance noted above for the solution f(z) of 



(6.1 ) under rescaling of the elliptic modulus k is indicative of a deeper invariance of 



20 



DMITRY KHAVINSON, ERIK LUNDBERG, AND RAZVAN TEODORESCU 



the solution: all the specific solutions in C discussed here are associated with fixed 
points in the moduli space of Riemann surfaces. 

Let again f(h(z)) be the analytic completion of a solution, and denote by Q the 
group of transformations which leaves supp(/x) invariant up to a global rotation. It 
follows that / is an automorphism of the quotient of the group of linear fractional 
transformations by Q, which can be in general a Kleinian group [30 . The limit set 
(accumulation points of the orbits of the group) can be finite (in which case it can 
consist of only 0, 1, or 2 points), or infinite. It is known ([2], Thm. 10.3.4.) that 
the set of homeomorphic solutions for a quasilinear elliptic equation of Laplace- 
Bcltrami type forms a group only in the case of finite limit set [2J. The Kleinian 
groups are called degenerate in this case, and they correspond to cither finite groups 
(with empty limit set), or the cyclic groups (generated by one element, with limit 
set consisting of 1 or 2 points). These correspond to the solutions described in the 
present paper (isolated point at infinity, respectively simple and double boundary 
point at infinity). 
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